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The  Finite  Sturm- Li ouvi lie  Transform 


By 

A.  Cemal  Eringen 
Illinois  Institute  of  Technology 

1.  Introduction 

It  is  known  that  certain  partial  differential  equations 
can  be  solved  with  the  use  of  particular  types  of  definite 
integrals  having  appropriate  kernels . The  choice  of  kernel 
depends  on  the  type  of  boundary  value  problem.  The  solution 
obtained  by  these  transforms  is  direct  in  the  sense  that  it 
contains  the  boundary  values  in  the  solution.  This,  of  course, 
is  lacking  in  the  classical  approach.  Finite  Fourier  trans- 
forms are  of  this  type  [1], 

Recently  Tranter  [2],  used  a Legendre  polynomial  as  a 
kernel.  Scott  [3]  following  a similar  approach  used  a Jacobi 
polynomial  as  a kernel  which  extends  the  result  of  [2], 

It  is  the  purpose  of  the  present  paper  to  extend  and 
unify  all  such  special  transforms.  Thus  we  employ  a kernel 
which  may  be  determined  to  suit  each  particular  type  of  problem. 
The  transform  can  be  employed  to  solve  a wide  class  of  linear 
second  order  partial  differential  equations. 

In  this  paper  we  deal  only  with  those  transforms  whose 
intervals  are  finite.  Thus,  the  results  obtained  here  are 
particularly  useful  for  finite  domains,  or  domains  which  are 
finite  in  one  direction.  An  extension  to  infinite  domains 
and  singular  cases  will  be  made  in  a later  paper. 
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2 


2 . The  Sturm- Li ouville  Expansion 
A second  order , linear,  homogeneous  differential 
equation  containing  an  arbitrary  parameter  \ has  the 
general  form: 

Mt  * poT,yy  * plT.y  * (p2  + Xp8>T  “ 0 (1) 

«0*y^bo 

where  ■ p^  (y)  5 v and  y are  the  dependent  and  real  inde- 
pendent variables,  respectively,,  Indices  after  comma  repre- 
sent differentiation,  i.e„:  4>  * d^/dy. 

a y 

Equation  (1)  can  be  transformed  into  the  canonical 


form: 


Lu  * \,u  , 


L s q(x)  --^-5  , a^z^b 
dx 


(2) 


where ; 

u - ©(x)v  , x -/( P3/pQ)1^2dy 

a ■ / 0(P3/P0)1,,2dy  , b - / °(p3/p0)1^2dy 

q(x)  - 6^f-  - *(x)  , 6(x)  - (p3k2/p0)1/4 

^(x)  - p2/P3  , k - exR/'(p1/p0)dy 


(3) 


It  is  simpler  tp  work  with  (2) . 

Ve  consider  two  solutions  ^(x,\)  and  ’^(x,*.)  of  (2) 
satisfying  the  boundary  conditions; 

l(a,\)  - sina  , 4 ,(a,\)  - -cos  a 

m A 

(4 

T^(b,\)  ■ sinp  , ^x(b,\)  - -cos  & 

Both  solutions  are  unique  if  q(x)  i£  real,  continuous  every- 
where in  (a,b)  and  has  finite  limits  for  x * a and  x ■ b 
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(see  [4],  p.6).  We  can  build  up  ^(x,\)  and  3^(x,>J  as  a 
linear  combination  of  two  independent  solutions  4>Q{x,\) 
and  V*  ,\)  of  (2).  For  let: 

Wo()J4(x,)J  - $0(x,\)[\0(aA)cos  a+  ^^x(aA)sina3 

- "^(x^C^U^cos  a+  ^x(a,>Jsina] 

<0 0(\>\<x,\>  - *0UA>[\0(bA)coSp+^0iJC(b,%.)sine]  ),  <5) 

- ’{Q(x,\)U>Q(b,\)cos  & + *0^x(b.)Jsinp] 

„0W  - W(*0>V  - *0lx'OWU''0  ‘ *o.xlx'Xl,Vx'»J  J 

Functions  £(x,\)  and  defined  by  (5)  satisfy  (2)  and  (4) 

Here  W(ufv)  is  the  Wronskian,  and  can  be  shown  to  be  indepen- 


dent of  x. 

Let  \n  be  a root  of  W(^/^)  ■ 0.  Then  we  have: 

“<V  - W‘  W *nW  *n.x\,  - 0 


or 


where: 

*„  * *(X>V  , and  \n  * 'Vx'xn>  • 

Hence  the  integration  gives: 

^(x#>.n)  - kn^(x,\n)  (6) 

Consequently  ^(x,^)  and  kn^(xf?^)  satisfy  all  of  the  boundary 
conditions  (4) . 

It  is  also  easy  to  show  that  ^(x,\n)  ijs  an  orthogonal 
set.  For,  in  view  of  (2),  (4)  and  (6)  we  have: 
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'VvAa.^  - • w(V*h)  - 0 

& a 

<w 

Hence  we  can  writes 


where  8^  is  the  Kroneclcer  delta  (equal  to  zero  f or  mfn, 
one  for  m*n) 

It  can  be  shown  also  that  all  roots  of  wtTv.)  * 0 
are  real  and  distinct  (see  [4]„  pp  11-12) . 

From  (5)  and  (6)  it  follows  thats 
tf0(bAn)cosp  ♦ ^o„x*be\i*sinP 
kn  * 0</a,\i^COSa  + Cx(a,Vs^na 

_ Vb-V costs  ♦ Ho.xto.^l^ne  (8) 

J^U.^looso  ♦ to  [xla>Vsino 

If  now  f (x)  is  an  inteqrable  function  over  (a,b)  and  if  a£x*b, 
we  have  the  Sturm-Liouville  expansion  of  f (x ) s 

f(x)  - E [k  /to  - (\  )]£(x,\  )/  0(y„\  ) f (y)  dy  (9) 

n.Q  n n n a n 

For  the  convergence  of  (9)  to  f(x)  see  [4]  pp  12-15  or  [5] 
pp275-276. 


3 . Finite  Sturm-Liouville  Transform 

Let  f (x)  be  a real  continuous  and  inteqrable  function  of  x, 
in  the  interval  (afb) „ We  define  the  finite  Sturm- Li ouv i 1 le 
transform  • Tta^)  associatec*  with  system  (2)  and  (4)  bys 

?(\n)  :3{f]-  / f (yltly.T^tdy 


(10) 


0 
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where  : 

t(y„V  - cy  vVl1/2*(r-V  (11) 

The  inversion  theorem  now  follows  from  (9) ; 

f(x)  -?T(V)t(x,v)  (12) 

n-o  n n 

Here  functions  ♦(ytf\n)  are  orthonormal,  i.eB; 
b 

S ^(x,*.  )t(x,\.  )dx  - S (13) 

a n a mn 

This  is  readily  seen  from  (7)  and  (11) „ 

Equations  (10)  and  (12)  are  basic  for  the  solution  of  some 
partial  differential  equations „ 


4.  Transform  of  Lf 


We  are  now  going  to  prove  thats 

SxlLf]  2 2x^q(x)f  • d2f/ax2}-  Bf(X.n,t) 

♦ \nT(\_,t)  (14) 

n n 


where : 


Bf<V*>  * Ckn/<oA(V]1/2  {-  flb.tjSfSi 


} 


(15) 


is  a function  depending  on  the  boundary  values  of  f(x,t), 
f (x,t)  at  x * a,b.  Subscript  f implies  that  the  function  B 

I* 

contains  boundary  values  of  the  function  f.  Here  Sx  repre- 
sents the  finite  Sturm- Li ovi lie  transform  taken  with  respect 
to  x: 
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To  prove  (14)  we  mutliply  Lf  by  t(x,\n)  and  integrate 
between  (a,b) , that  is,  by  definition  we  haves 

S^Lfj  * [q(x)f(x,t)  - a2f/ax2]t(x,\n)dx 

Integrating  the  second  term  in  the  integrand  of  (17)  twice 
by  part  we  obtain: 

S^Lf  } - -[f^x  t - ^xf]^  + f(x,t)Lt(xAn)dx 

If  we  now  use  (4)  and  (6)  in  the  first  term  of  the  right 
hand  side  and  (2)  inside  the  integrand  we  obtain  (14) . 

5.  Solution  of  Some  Partial  Differential  Equations 

a)  Let  us  consider  a parabolic  partial  differential 
equation 

Mv  ♦ p^v  i • 0 

where  Mv  is  defined  by  (1),  p^  * (y#t)  are  known  and 

v(y,t)  is  the  unknown  dependent  variable. 

The  use  of  transformation  (3)  reduces  (19)  to: 

u,xx  * [*•  " P(*.t)]u  + (p4/p3)u^t  " 0 

where: 

P (x, t ) - q(x,t)  - p4p0(  k/2po)2  ^|>(p0/p3  k2) 

Here  q(x,t),  k(t)  and  \{t)  are  given  by  (3)  with  t as  a 
parameter.  If  $(x,t,>.n)  and  t(x,t,\R)  are  now  obtained  as 
before  except  with  q(x)  replaced  by  p(x,t),  t a parameter. 


(16) 


(17) 


(18) 


(19) 


(20) 


(21) 


- 7 - 

and  P4/P3  a function  of  t alone,  then  (20),  with  the  use 
of  the  finite  Sturm-Liouville  transform,  can  be  reduced 
to  an  ordinary  differential  equation  in  t.  For,  let 
P4/P3  ■ r(t),  then  application  of  the  transform  to  (20) 
gj.ves: 

-Bu(Vt)  + (**“>.n)u(>.n,t)  ♦ r(t)u  “ 0 (22) 

Equation  (22)  is  of  first  order B The  complete  solution  is; 

u(Vt)  - C(\n)exp(/t  n “ dt)  + [exp  (/*— n dt)]  » 

/t(Bu/r)  [exp  X dt)]dt  (23) 

Now,  inversion  theorem  (12)  gives: 

u ( x , t ) - E u(\n,t)tn(x,t,\n)  (24) 

n-o 

Therefore,  the  solution  of  (19)  is  effected,,  We  must,  however, 
remember  that  at  least  one  of  each  pair  of  boundary  values 
u(a,t),  u (a,t)  and  u(b,t),  u (b,t)  must  be  known.  By 

9 * 9 * 

selecting  a and  3 properly  in  Bu(*-n,t)  we  can  make  the  terms 
containing  the  unknown  pair  zero.  Of  course  any  knowledge  of 
the  boundary  values  of  u leading  to  the  complete  evaluation 
of  Bu  is  sufficient  for  the  solution. 

b)  A second  order  linear  partial  differential  equation 
containing  a parameter  \ has  the  general  form: 

Nv  ♦ a4v  z ♦ a&v  T ♦ (ag  ♦ \a7)v  • 0 
Nv  • a,v  ♦ 2a«v  ♦ a«v  ^ 

i , ZZ  l ,ZT  <3  , TT 


(25) 
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where  a^  * a^(z,x).  If  the  equation  is  of  elliptic  type 
we  have: 

ala3  " a22  ^ 0 (26) 

Consider  now  an  elliptic  equation  (25).  Let  £(z,x)  * 
const,  and  t)(z,t)  * const 0 „ respectively  be  the  solutions 
of  the  following  differential  equations: 

dx/dz  - C^z^t)  dx/dz  - C2(z„x) 

**l\  a2  + i[— - 

c2J  ~*im  &1  ai 

If  we  set: 

2y  “ l + *)  , 2it  - l - t)  (28) 


equation  (25)  can  be  transformed  into  the  canonical  form  [6]: 
\yy  * pl\y  * (p2  * **p3)v  * p4\t  * \tt  ■ 0 (29) 


where: 


PX  * f py  » P2  “ 2a6^  ’ p3  " 2a7^ 


fpt  . S>  - a^ (y 


J * t h + 2a«  (y  y , + t .t  „) 
,Z  , Z £ »z  ,x  *.  T 


♦ a3(y^2  + t'XZ)  (30) 

Here  operator  P is  defined  by: 

P*  ■ + a4^#z  + a5^fx  OD 

and  is  given  by  the  second  of  (25). 

Using  the  transformation  (3)  with  pQ  ■ 1,  (29)  becomes: 

U,XX  + <>•  - p)u  + rlu,t  + rou,tt  ' 


0 


(32) 


where : 
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p - q(x,t)  - mp4g  . - , rQ  - P3-1 

rx  - mp4g  + 2mg^t  , g - (p^2)"1/4' 

m-  (k2/^3)1'4’  (33) 

In  the  special  case  where  r^,  and  rQ  ^re  functions  of  t alone 
equation  (32)  can  be  transformed  into  an  ordinary  differential 
equation  by  applying  the  finite  Sturm- Li ouvi lie  transform. 
Hence: 


rou,tt  * rlu,t  * - Bu(\n,t)  (34) 

We  can  now  either  use  transformation  (3)  with  respect  to  t, 
with  xQ,  r^  replacing  pQ,  p^  and  p2  * 0,  p3  ■ 1 and  find 
the  solution  of  (34)  and  then  invert  it  or  we  can  solve  (34) 
directly  for  u(\n,t)  and  then  (24)  gives  u(x„t). 

It  must  be  remembered  that  this  method  has  the  advantage 
that  it  contains  the  boundary  values  of  the  function  within 
the  solution,  consequently,  the  difficulty  of  satisfying 
boundary  conditions  is  eliminated.  The  method  is  also  valid 
in  the  cases  where  the  usual  technique  of  separation  of 
variables  fail. 

6.  Solutions  and  of  Some  Special  Equations 

Below  are  given  two  independent  solutions  and  )^0 
of  some  special  differential  equations  (D.  E.)  of  mathe- 
matical physics.  Their  canonical  forms  (2)  will  be  given 
by  only  determining  u,  q,  and  X. 
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a)  Bgasel  Functions 

v + i y 4 (s2  - 4)v  - 0 (D.  E„) 

,yy  y \y  V8  u vu* 

u - x^v  , y - x , q “ (v2  - ^)x“2  \ - s2 

*0(x,\)  - x1/2  Jv  (xs)  , ^0(x,\)  ■ x^^fxs) 

Interval  (a,b)  should  not  contain  x * 0. 

b)  Spherical  Harmonics 

2 ..2 

C(l-y  )▼  J]  + [v(v4l)  - — 9]v  « 0 (D.  E.) 

#y  .y  lV 

2 

u * v cos1/^  y « sin  x , q - — - i tan2  x - i , 

cos  * ’ 

X • v (v  + 1)  f ^0(x,\)  * cos1/2  x • P^tsin  x ), 

)^0(x#\)  - cos1/2  x "Q^ ^ (sin  x ) 

Interval  (a#b)  should  not  contain  (-1,  1). 

c)  Her  mite  Polynomials 

\yy  - yv#y  + nv  - 0 (D.  E.) 

.-2ii  i —2 

u ■ e x ' 4 v#  y-x  , q«“7*T-  , * n , 

2 2 
$0(x,\)  * e"x  /*  Hen(x)  » ?^0(x,\)  “ e‘x/4hen(x) 

d)  Tschebyscheff  Polynomials 

(i-y2W  ^ “ yv  + n2v  - 0 (D.  E.) 

2 

u - v , y ■ sin  x,  q ■ 0 , \ ■ n 

#0(x#\)  - Tn(sin  x)  - cos  n(sin“1y)  , 

5^c(x,\)  ■ Un(sin  x)  ■ sin  n (sin*  "Sr) 
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e)  Mathieu  Functions 


u v_  ♦ (X.  - 2h^coe2x.(U  * 0 

f A3fc 


<D.  E.) 


V1-0  ’ °»2n(x)'  oe2n+l(x)  • se2n<x)>  s2n+l(x) 


V*'70  - ce2n(2)(x),  082n+1l2,(x),  ae^fx) 


(2) 


.(2), 


8 (x) 

s2n+l  w 


(h  * 0,  1,  2,  ...  except  for  se2n  only  n ■ 1,  2,  ..„) 
>0(x,\)  is  periodic. 

f)  Wittaker  Functions 


1 - n2 

V^yy  + (-^  + J + )V  " 0 <D*  E.) 

u » (x/2)^2v,  y » x2/4,  q * ♦ y*. 

4x^  D 

*QU,X.)  - (x/2)1/2  Wk^(x2/4),  )^(x,X.)  - (x/2) 

W (-x2/4)  (*) 

*k,H 


1/2 


Interval  (a,b)  should  not  contain  x * 0. 


7.  Application 

As  an  illustration  we  solve  a heat  conduction  problem 
which  may  be  looked  upon  as  a mathematical  model  of  volcanos. 
The  problem  also  has  application  in  the  exhaust  parts  of  jet 
engines.  As  far  as  I know  the  problem  has  not  previously 
been  solved. 


Here  comma  does  not  mean  differentiation 
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Consider  the  conical  sheel  enclosed  by  two  coaxial 
cones  having  the  same  apex  0(  and  two  concentric  spheres 
having  0 as  their  center  (Fig.  1) . 


In  polar  coordinates  the  inner  and  outer  conical 
surfaces  are  given  by  0 ■ 0Q  and  9 • 0^,  and  the  end 
surfaces  by  r ■ rQ  and  r • r^. 

Problem:  Determine  the  steady  temperature  distribution 

within  the.  shell  under  the  general  axi-symmetric  boundary 
conditions  in  temperature.  This  is  the  general  Dirichlet 
problem  for  the  domain  under  consideration.  The  differential 
equation  and  boundary  conditions  (B.  C.)  are  given  below: 
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(D.  E.) 


(B.  C.) 


AV  - r~2  (r2V  ) ♦ r~2[  (l-y2)V  1 v - 0 

ff1  9J  9 J 

y - cos©  ,rc<r<r1  , ©0<6 

V - V (r)  for  6 - 6^ 

o o 

v - ^(r)  for  e - e2 

V - V3 (©)  for  r - rQ  \ 

l eft<e<e1 

V - V4(6)  for  r - r1  J ° 1 


rft<r<ri 


(35) 


(36) 


Here  V(r„y)  is  the  temperature  function,  r and  0 are  the 
polar  coordinates. 

Solution:  We  can  exclude  the  second  term  of  (D.  E.)  (35) 

if  we  use  a finite  Sturm- Li ouvi lie  transform  associated  with 
the  Legendre  equation: 

[ (l-y2)v  ] + v(v+l)v  - 0 (37) 

9 y o y 


In  view  of  (6f  b)  with  jx  ■ 0 we  find  that  if  we  select 
y ■ sin  x,  u * v cos^2  x we  can  transform  (37)  to  canonical 
form  leading  to  solutions: 

$Q(x,\)  - cos^2x  " Py(sin  x)  , )^(xf\)  “ cos1^2  x 


° (sin  x)  y \ » v'v  + 1)  (38) 

where  Pv  and  are  Legendre  functions  of  the  first  and 
second  kind,  r espectivelya  In  view  of  (14),  when  we  apply  the 
transform  to  (35)  the  second  term  gives  - vn(vn+D  “ Bv(r,\n)« 
Now  By  contains  four  arbitrary  functions  which  must  be 
specified  on  the  surfaces  0 ■ 0Q  and  0 ■ 0^.  The  terms  con- 
taining the  derivatives  of  v are  not  given B Thus  if  we 


- 14  - 


l 

£ 


select  a • 8 ■ 0 these  terms  drop  out,  leaving  the  terms 
containing  v0(r)  and  v^(r),  which  are  given.  Hence 
^(x,\)  ando>0(X)  of  (5)  become: 

4(x,\)  - cos^2x0  cos^2x  [Pv (sin  xJQ^ (sin  xQ) 


- Py (sin  x0)Qv (sin  x)]  , 


- S - 6 


- e. 


o>0  <\) 


X1  7 " "1  * 

By  (11)  and  (8)  we  have: 

- Ckft/®^^)]1^2  *<*»**) 

,1/2x,  P (sin  x,) 

~ vn 


cos 


n 


<=os1/2x0  Pv  (sin  x0) 
n 


Calculation  of  <o(\)  m t ~ x gives 

00 (X.)  * cos^2xo  cos1^2x1  [Pv  (sin  xq)Qv  (sin  x-^) 
- Pv(sin  Xj^JQy  (sin  xQ)] 


(39) 


(40) 


(41) 


Consequently,  the  roots  • vn(vn+l)  of  ®*(\)  *0  satisfy 
the  following  equation 

Pv  (sin  Xj^)  Qy  (sin  x^ 

Pv isin  x0)  “ Qy  (sin  xQ)  (42 

Therefore  ♦(x,*^)  is  completely  determined. 

After  transforming  (35)  with  y ■ sin  x and  u “▼  cos^2x 
we  apply  the  finite  Sturm-Liouville  transform.  That  is,  we 
multiply  the  equation  by  tU,*^)  and  integrate  between  xQ  and 
The  result  is: 
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* W1**-  Vr*V 

Vr*V  •tknfcAVl1,J[V„( r)  - k^d)] 

Equation  (43)  is  an  ordinary  differential  equation  of  Euler 
type  whose  solution  can  be  found  by  variation  of  parameters. 
Hence: 

-v  -1 


tf(rfvn)  -C1(vn)r  n ♦ C2(vn)r  n + F(r„vn) 

_ Bv  (p  ) v -v  -1,  -v  -1  V , 
F(r„v  ) mS  —y+'kr-  <r  P -r  n p 


)dp 


(43) 


(44) 


Let  the  transforms  of  V«  (0)  and  V.  (0)  be  V«(v„)  and  V.(v  ). 

3 4 3 n 4 n 

Cj  and  C2  will  then  be  determined  from  the  remaining  conditions: 


V - V3(vn>  , for  r * r( 


e. 


„ for  r - r^ 


,<e<e1 


This  gives  two  linear  equations  for  Cj  and  C2  whose  solutions 

are: 


[^(vn)-P(>0.vn)]r7n~1-L^(vn)-r(r1.vn)]r0'Vlt~1 

Cl(vn)-  r-v„-l  ~ “ 

ro 


’“-T 

■X1  n r.  n 


(45) 


- ...  CW-F(rl'VroVn  -t^(vn)-P(r0.vn)]riV» 
2 n **  =v* 

r n r,  n -r,  n r n 


r n r, 

Oi  JL  U 

Hence,  V(r,vR)  is  completely  determined.  The  inversion  theorem 
(12)  now  gives: 

V(r,x)  - l f (r,*B>*(*.)Sl>  (46) 

where  the  summation  is  extended  over  all  roots  of  (42) . 


- 16 


The  analysis  given  above  is  formal.  It  can,  however, 
be  made  rigorous  by  showing  that  the  solution  satisfies  both 
the  differential  equation  and  the  boundary  conditions. 

It  may  be  worth  while  to  remark  that  with  the  use  of 
the  method  of  separation  of  variables  the  solution  of  the 
above  problem  would  have  been  difficult  and  lengthy,  if  not 
impossible. 
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